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ABSTRACT 


.  A  model  for  the  deformation  of  multi-layer  metallic  laminate  targets 
during  ballistic  impact  is  described.  The  model  divides  the  process  into  a 
stage  of  plug  acceleration  followed  by  dishing  of  the  rear  of  the  target,  and 
allows  for  the  effects  of  variable  interlamellar  bend  strength.  A 
computational  procedure  is  described,  as  well  as  the  mechanical  test 
procedures  to  generate  the  input  data  for  the  program.  Examples  of  the 
application  of  the  model  are  described  to  illustrate  the  interpretation  of 
results,  the  limitations  of  the  model,  and  the  use  of  the  model  is 
parametric  studies  for  laminate  design.  The  computer  program  is  listed 
with  typical  input  and  output  data.  f 
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A  COMPUTATIONAL  MODEL  OF  THE  PERFORATION  OF 
MULTI-LAYER  METALLIC  LAMINATES 


1.  INTRODUCTION 

Computation  of  ballistic  impact  events  can  be  via  simple  analytical  equations, 
more  complex  analytical  models  requiring  computational  procedures,  or  by 
finite  element  or  finite  difference  computer  codes.  Whilst  the  computer  codes 
have  been  partially  successful  in  providing  an  impressive  graphic  output  of 
stages  in  the  perforation  process,  they  have  not  been  able  to  effectively 
provide  a  means  of  calculating  the  result  of  ballistic  impact  events.  The 
analytical  equations  and  procedures  have  provided  the  main  tool  for  efficient 
parametric  studies  and  basic  insight  into  the  mechanics  of  perforation,  provided 
the  scope  of  the  problem  was  well  defined.  A  number  of  such  analytical 
procedures  have  been  developed  and  documented  in  MRL  work  11-4),  and  it  has 
been  shown  that  an  understanding  of  the  mechanics  allows  the  methods  to 
provide  good  approximate  answers  to  a  range  of  problems  [41. 

In  recent  work  Woodward  and  Crouch  14]  demonstrated  how 
approximate  calculations  could  be  performed  on  simple  two  layer  laminates  by 
inserting  small  modifications  in  an  earlier  plugging  failure  model  [11.  In  the 
present  work  the  same  general  approach  is  used  with  experimental  observations 
on  the  behaviour  of  multi-layer  metallic  laminates  to  develop  a  model  for 
perforation  of  such  laminates.  The  step  from  homogeneous  metals  to 
laminates  also  required  the  development  of  suitable  test  procedures  designed  to 
provide  appropriate  materials  input  data. 


2.  PERFORATION  MODEL 

For  practical  purposes  the  number  of  layers  in  the  metallic  laminates  is  defined 
as  the  number  of  layers  of  metal,  so  that  any  adhesive  or  fibre  reinforcement  is 
considered  as  part  of  this  unit.  The  thickness  of  a  layer  is  then  the  total 


thickness  divided  by  the  number  of  metal  layers,  and  the  approach  assumes  all 
layers  are  equivalent.  The  problem  to  be  considered  assumes  a  non-deforming 
projectile. 


The  perforation  process  is  divided  into  two  stages  as  illustrated 
schematically  in  Fig.  1.  When  the  projectile  impacts  the  target  the  velocity 
difference  ensures  an  initial  stage  of  acceleration  of  material  in  front  of  the 
projectile.  During  the  time  in  which  the  plug  is  accelerated,  the  reaction 
force  slows  the  projectile  so  that  a  stage  is  reached  after  which  projectile  and 
plug  move  at  the  same  velocity.  In  a  homogeneous  metal  target,  shearing 
between  the  plug  and  the  target  would  then  lead  to  fracture  and  separation  of 
the  plug.  In  a  laminate,  on  the  other  hand,  the  difficulty  of  crack  propagation 
across  a  number  of  laminae,  combined  with  the  ease  of  separation  by  a  crack 
travelling  along  the  bond  ensures  that  the  second  stage  of  perforation  involves 
delamination  and  dishing  failure  of  the  rear  of  the  target  as  illustrated  in  Fig. 
1(c).  This  second  stage  must  involve  an  integral  number  of  layers  as  it  is 
assumed  that  the  bonding  between  lamina  is  weak  compared  to  the  metal 
strength. 


(a)  Stage  I  -  Plug  Acceleration 

The  method  follows  generally  the  procedure  described  by  Woodward  and 
de  Morton  111  with  modifications  allowing  adaption  to  composite  materials  for 
which  flow  rules  have  not  been  universally  established.  The  uniaxial  stress  and 
uniaxial  strain  data  from  a  constrained  compression  test  on  the  composite, 
described  in  Section  3,  are  curve  fit  to  the  form 


^2 

o  =  o~6  ,  e  >  e  (lb) 

2  o 

where  a  is  the  normal  stress  on  a  punch  in  the  constrained  compression  test, 
e  is  the  uniaxial  strain, 

o  ,  o  ,  rj  ,  n„  and  e  are  curve  fitting  constants  for  data. 

The  use  of  two  curves  separated  at  a  strain  £q  allows  considerable  flexibility  in 
covering  a  wide  range  of  material  behaviour.  The  elastic  phase  of  loading  is 
not  explicitly  separated  as  this  leads  to  an  insignificant  change  in  the  result. 
The  assumed  uniform  velocity  of  a  plug  moving  ahead  of  the  projectile  is 
obtained  by  integration  of  the  wave  equation 


AV 


A  o 
pc 


(2) 
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where  Av  is  a  velocity  increment  associated  with  a  stress  gradient  Ao  and 
p  and  c  are  the  material  density  and  wave  velocity  respectively. 


This  plug  velocity,  V1(  becomes 


v  _  (Vl,l/2  __2 _  ^Hj  +  D/2 

i  p  n,  +  i  e 


,  e  &  e 


(3a) 


and 


Vl  1/2  _ 2  (’)1+1^2  -'72°2.1/2  2  (n2+1)/2  (')2+1)/2 


V1  (  P  ’  Hj+X  60 


+  (-*-*)  . 
p  n  o+l 


(e 


) ,  e  >  e 


(3b) 


The  normal  force  F  on  the  plug  of  diameter  D  is  given  by 


F  =  \  D2  oxz  2  ,  e  s 

F  =  oJ2 


4  2 


o„e  ,  e  >  e 


(4a) 

(4b) 


When  the  back  of  the  plug  moves  a  distance  ALg  an  increment  of  work  aW^  is 
done  where 


-  fai^ 


(5) 


At  the  strain  e ,  the  plug  height  is 


h  =  hQ  exp  ( -  e ) 


(6) 


where  hQ  is  the  initial  laminate  thickness. 

The  situation  at  the  end  of  stage  I  is  illustrated  in  Fig.  2(a)  which 
relates  the  displacement  of  the  back  of  the  plug,  Lg,  the  front  of  the  plug,  Lp, 
the  residual  contact  length,  Lq,  and  the  initial,  hQ,  and  final,  h,  plug 
thicknesses  to  both  the  geometry  used  in  the  plugging  solution  of  Woodward  and 
de  Morton  [11  and  the  shearing  of  the  laminae  in  the  present  instant.  The 
following  relations  follow  directly  from  the  earlier  plugging  model  til. 


t  -  2(h0-h)/V0 


(7) 


where  t  is  the  time  for  stage  I  and 

VQ  is  the  projectile  impact  velocity, 

Lp  =  Vj  t/2  (8) 


Lg  =  (VQ  +  V1)t/2 


(9) 


The  work  done  in  shear,  Wg,  in  stage  I  is  given  by 


Wg  —  ttD  {h(hQ  —  h)  +  L^,(2h  —  Lp)}  (10) 

where  t  is  the  shear  strength  of  the  laminate.  The  work  done  in  friction 
between  the  projectile  and  the  laminate,  Wp,  is 


WF  =  ttD  t  Lb2/2 


(11) 


The  new  projectile  velocity,  Vp,  at  the  end  of  stage  I  is  given  by 
Vp  =  (VQ2  -  2  (WN  +  Wg  +  WF)/m>1/2 


(12) 


where  m  is  the  projectile  mass. 

The  computational  procedure  increments  strain  in  the  plug  in 
stages,  calculating  new  values  of  (plug  velocity)  and  Vp  (projectile  velocity) 
until  a  value  is  found  at  which  the  work  done  is  sufficient  for  Vp  to  equal  Vp 
At  this  stage  we  have  a  situation  where  projectile  and  plug  are  moving  together 
and  we  proceed  to  stage  n. 


(b)  Stage  II  -  Dishing  Failure 

For  a  homogeneous  plate  where  plugging  failure  is  observed,  Woodward  and  de 
Morton  (11  used  a  shearing  failure  mode  for  exit  of  the  plug.  For  the 
laminates,  some  mechanism  which  allows  dishing  and  stretching  is  desirable, 
and  for  the  particular  case  of  flat  ended  projectiles  plug  separation  as  well. 
Two  decisions  need  to  be  made,  viz.  how  many  laminae  remain  to  be  perforated 
after  stage  I,  and  which  model  to  use  to  calculate  the  work  done. 

Examination  of  a  number  of  sectioned  partly  perforated  samples 
was  used  to  try  and  establish  a  criterion  for  which  laminae  were  sheared 
through  to  failure,  and  which  remained  intact.  It  was  generally  observed  that 
significant  shear  displacements,  greater  than  a  layer  thickness,  are  possible 
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without  separation  occurring.  In  addition  shear  zone  widths  vary  from 
material  to  material  with  the  distinction  between  shearing  and  bending 
deformation  becoming  one  of  judgement  in  many  cases.  Thus  a  simple 
displacement  or  shear  strain  criterion  for  failure  was  not  seen  to  be  possible  at 
this  stage.  Based  on  observations  it  was  decided  to  pragmatically  locate  the 
back  of  the  plug  at  the  end  of  stage  I,  defined  by  Lg  and  the  next  interface 
below  this  point  as  the  one  at  which  separation  occurred,  as  illustrated  in 
Fig.  2(a)  and  (b).  This  approach  concurred  well  with  what  was  observed  in  most 
accompanying  experiments  and  the  ballistic  limit  predictions  were  reasonable. 

For  a  model  of  deformation  in  stage  II  the  structural  approach  of 
Woodward  [31  appears  ideal  as  it  allows  dishing,  stretching  and  plugging 
deformations.  Unfortunately  the  assumptions  of  that  model  apply  only  in  the 
case  of  thin  targets  where  projectile  radius  is  greater  than  the  target  thickness, 
and  when  the  target  thickness  increases  beyond  this  the  model  produces  large 
errors.  For  many  problems  of  the  type  being  considered  in  the  present  instant 
it  was  observed  that  after  the  first  stage,  the  thickness  still  to  be  penetrated 
was  generally  of  the  order  of  one  projectile  diameter,  thus  the  use  of  the 
structural  model  is  not  possible.  There  are  no  other  models  which  treat  all 
aspects  of  the  deformation  and  which  at  the  same  time  do  not  require  ballistic 
experiments  to  generate  some  of  the  input  deformation.  For  many  blunt 
projectile  impact  problems  it  has  been  found  [31  that  a  simple  analytical 
equation  for  dishing  failure  which  allows  for  bending  and  stretching  failure  only 
of  the  plate  gives  good  estimates  of  the  work  done.  This  model  was  originally 
derived  by  Taylor  [51  and  Thomson  [61  and  modified  by  Woodward  [71  to  include 
the  effect  of  plate  bending.  The  work,  WD,  is  given  in  terms  of  the  thickness 
Te  to  be  penetrated  in  stage  n.  Fig.  2(a),  and  the  material  strength,  Y,  by 


W  =  n DYTg  (D  +  "Te)/8 


(13) 


where  the  first  term  represents  the  stretching  work  and  the  second  term  the 
work  done  in  bending. 

For  want  of  a  better  approach  this  method  of  treating  the  dishing  failure  was 
adopted  as  it  includes  some  of  the  observed  features  and  should  calculate  the 
work  done  to  a  reasonable  approximation.  Whilst  it  copes  with  these  aspects  it 
must  be  accepted  that  it  may  not  be  able  to  distinguish  differences  in 
performance  which  result  from  small  changes  in  failure  mode.  To  cope  with 
the  differences  in  bend  properties  of  a  laminate  and  a  homogeneous  metal  and 
account  somewhat  for  interlamellar  bond  strength  effects,  equation  (13)  was 
modified  to  accept  appropriate  laminate  strength  terms. 

If  there  are  n  layers  to  be  perforated  and  the  layer  thickness  is  TE, 
then  the  tension  or  stretching  work,  Wj,  in  the  second  stage  is 

WT  =  n  D2  o  nT.  IS  (14) 

*  u  L 
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where  a  is  the  ultimate  tensile  strength  of  the  laminate  measured  in  a  tension 
test. 

In  bending  the  problem  is  likened  to  a  beam  of  width  77D  bent 
through  an  angle 

77/2.  The  beam  is  assumed  to  be  fully  plastic  and  0 p  is  the  fibre  stress 
measured  in  a  simple  four  point  bend  test.  Then  for  an  unbonded  composite 
the  work  done  in  bending  is 


WB  =  nrr 


0fOTl2/8 


(15) 


For  a  bonded  laminate  bending  is  divided  into  stages  before  and  after  the  bond 
breaks  along  the  centreline  at  some  angle  of  bend  e  given  by 


e 


e 

c 


(h  /nT.  )2 
O  L 


(16) 


where  e  is  the  angle  at  which  the  bond  breaks  in  .he  simple  four  point  bend 
test  usecl  in  the  original  laminate  of  thickness  hQ  to  determine  Op. 

The  work  done  in  bending  for  a  bonded  laminate  is  then  given  by  the 
sum  of  the  two  terms 


WB  =  n2  rr  ap  D  T^e/4 
and  WB  =  (nx2  +  n22)  injpD  (|  -  <?)/ 4 
where  the  integers  n^  and  n2  are  defined  by 


(17a) 

(17b) 


nl  +  n2  =  n 


(18a) 


and  I n^  -  n j !  si 


(18b) 


which  allows  for  the  n  to  be  uneven  so  that  the  laminate  effectively  breaks  its 
adhesive  bond  along  the  nearest  bondline  to  the  centre. 

The  solution  procedure  is  to  take  the  kinetic  energy  of  tne 
projectile  and  plug  from  stage  I,  calculate  the  work  done  in  stage  II  using 
equations  (14)  to  (18)  and  compare  these  to  determine  whether  or  not  the 
projectile  perforates.  The  computer  program  to  carry  out  the  calculations  was 
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written  so  that  it  increments  the  initial  velocity  in  small  steps  until  the 
projectile  just  perforates  and  this  is  treated  as  a  prediction  of  the  ballistic 
limit  of  the  composite. 


3.  LAMINATE  CHARACTERIZATION 

Use  of  the  above  model  requires  materials  data  in  a  suitable  form  to  treat  plug 
compression  and  acceleration,  shear  strength,  bending  and  tensile 
deformation.  With  the  development  of  a  better  understanding  of  composite 
behaviour  it  may  be  possible  to  estimate  these  parameters  from  the  properties 
of  the  component  parts.  For  the  present,  a  set  of  procedures  has  been 
developed  which  ‘simulate’  the  deformations  experienced  in  the  ballistic  impact 
as  best  possible,  so  that  the  data  is  appropriate,  whilst  the  testing  remains 
simple.  Many  of  the  parameters  are  combinations  of  basic  strength  data,  such 
that  experience  with  application  to  a  number  of  composites  may  allow 
estimates  of  the  parameters  with  far  more  limited  testing  programs  in  the 
future.  Rate  dependence  of  mechanical  properties  is  not  considered  as  the 
aluminium  alloys  used  in  the  present  work  are  insensitive  to  strain  rate.  There 
is  however,  no  difficulty  in  including  tha  effect  of  strain  rate  on  strength 
properties  in  the  analysis,  if  that  data  are  available  The  procedures  are 
described  below. 


(a)  Compressive  Deformation 

The  compression  of  the  plug  during  acceleration  can  be  simulated  by 
compression  of  a  plate  between  two  flat  axially  aligned  punches,  a  constrained 
compression  test.  The  subpress  for  carrying  out  such  tests  was  described  by 
Woodward  and  de  Morton  [11,  and  its  application  to  laminates  illustrated  by 
Crouch  and  Woodward  [81.  Typical  uniaxial  stress/uniaxial  strain  data  from 
such  a  test  is  shown  in  Fig.  3(a)  with  a  typical  cross  sectioned  sample  in  Fig. 
3(b).  The  data  is  easily  fitted  to  curves  of  the  form  of  equations  (1)  to  obtain 
the  materials  parameters  o^,  rj^,  a r /  and  e  ,  using  straight  lines  on 
plots  of  logarithm  of  stress  versus  logarithm  of  strain. 


(b)  Shear  Deformation 

For  shear  deformation  the  laminate  is  assumed  to  approximate  a  rigid/plastic 
material  with  constant  flow  stress.  A  blanking  test  is  carried  out  to  produce  a 
force  displacement  graph  as  in  Fig.  4.  The  area  under  the  curve  is  integrated 
up  to  a  displacement  equivalent  to  the  laminate  thickness.  From  the  work 
done  in  shear,  Wg,  the  punch  diameter  D,  and  the  laminate  thickness  hQ,  the 
she  ir  stress  for  use  in  the  program  r  ,  becomes 

r  =  2  W„/t7D  h  2  (19) 

s  o  o 
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This  is  equivalent  to  assuming  a  load/displacement  relation  of  constant  slope 
with  total  work  up  to  the  displacement  of  the  sample  thickness  equal  to  the 
actual  work  in  the  test  as  shown  in  Fig.  4. 

(c)  Bending  Deformation 

Figure  5  shows  a  set  of  typical  moment/angle  of  bend  curves  obtained  in  four 
point  bending.  In  the  case  of  well  bonded  laminates,  the  lamina  remained 
bonded  to  a  very  large  bend  angle.  With  "weak"  bonding  a  discontinuity  is 
observed  in  the  moment/angle  curves  when  bond  failure  occurs  along  the 
centreline  of  the  laminate.  Generally  no  further  bond  breakage  was  observed 
to  very  large  bend  angles,  so  beyond  the  discontinuity  the  sample  can  be 
considered  as  being  made  up  of  two  pieces  of  multilayer  laminate  with  no  bond 
in  the  interface  between  the  two.  For  unbonded  laminates  the  bending 
moments  were  consistent  with  the  bending  of  the  sum  of  the  individual  layers 
of  the  laminate.  Thus  from  the  moment  to  which  the  curves  extrapolate  at 
high  angles  of  bend,  Mp,  the  fibre  stress,  Op,  was  determined  from  equations 
(20a)  for  good  bonding,  (20b)  for  "weak"  bonding  and  (20c)  for  unbonded 
laminates. 


=  4  M  /w  h 

P  o 


(20a) 


8  M  /W  w  2 
P  ° 

4  n  M  /w  h  2 
P  o 


(20b) 

(20c) 


where  w  is  the  width  of  the  sample  used  in  the  four  point  bend  test. 

The  critical  angle  of  bend  at  which  the  bond  breaks  in  the  four  point 
bend  test,  e  ,  was  also  measured  and  used  as  a  measure  of  the  strength  of  the 
bond  in  the  computations  for  the  cases  of  "weak’’  bonding. 


<d)  Tensile  Deformation 

The  parameter  used  for  the  tensile  flow  stress,  o  ,  in  equation  (14)  was  the 
ultimate  tensile  strength  measured  in  a  tension  test  on  the  composite. 


(c)  Experimental  Laminates 

The  construction  and  physical  properties  of  laminates  tested  in  the  present 
study  are  given  in  Table  I,  and  the  input  data  for  simulation,  determined  with 
the  above  testing  procedures  is  given  in  Table  II . 
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4.  DISCUSSION 


Many  of  the  assumptions  in  the  simplified  model  have  been  explicitly  stated 
above.  Three  assumptions  are  based  purely  on  practical  approaches  to 
observed  behaviour,  viz.  the  method  used  to  estimate  the  number  of  layers  still 
to  be  perforated  at  the  end  of  stage  I,  the  simplified  method  of  including  bond 
strength  as  measured  in  the  bend  test,  and  the  use  of  a  plug  mass  in  stage  n 
rather  than  the  mass  of  a  whole  dished  region  as  the  latter  is  rather 
indeterminate.  The  treatment  of  stage  n  in  terms  of  the  simple  analytical 
dishing  model  gives  reasonable  estimates  of  work  done,  however  it  will  prevent 
the  model  distinguishing  material  differences  which  relate  to  the  ductility  of 
the  plate  as  seen  in  the  case  of  plug  separation.  It  has  been  observed  that  the 
model  predicts  lower  than  expected  ballistic  limits  with  unbonded  laminates 
when  impacted  by  flat  ended  projectiles  relative  to  bonded  laminates.  The 
difference  in  experiments  is  not  generally  that  great  and  the  discrepancy  arises 
because  the  model  attributes  more  work  to  the  bending  deformation  of  the 
bonded  compared  to  the  unbonded  layers.  In  practice  an  unbonded  laminate 
can  show  greater  ductility  in  bending  and  this  aspect  is  not  part  of  the 
formulation. 

What  then  is  the  application  of  the  model?  Table  IH  presents 
observed  and  predicted  ballistic  limits  for  a  range  of  laminates  impacted  by 
flat  ended  projectiles.  The  cases  include  variations  in  bond  strength,  the 
inclusion  of  Kevlar  fibre  reinforcement,  different  layer  numbers  and 
thicknesses  of  laminate  as  well  as  unbonded  laminates.  Generally  the  estimate 
of  the  ballistic  limit  is  of  the  correct  order  of  magnitude.  In  addition  where 
conditions  are  similar  the  order  of  merit  is  correctly  predicted  in  most  cases. 
The  model  is  seen  to  generally  overestimate  performance.  Cases  D1  and  Ml 
illustrate  the  case  of  an  unbonded  laminate  (Dl)  showing  higher  ballistic  limit 
relative  to  the  bonded  laminate  than  expected  because  of  greater  ductility, 
despite  the  expected  much  higher  ballistic  limit  in  the  bonded  laminate  (Ml) 
which  has  a  higher  bond  strength.  The  low  ballistic  limit  in  the  case  of  M3  is 
associated  with  poor  ductility  leading  to  complete  fracture  separation  of  petals 
at  the  back  of  the  plate.  Figure  6  compares  a  section  from  a  perforated 
laminate  made  from  the  2090  alloy  with  a  more  ductile  failure  from  a  laminate 
of  2024  alloy  and  a  typical  six  ply  7075  aluminium  laminate. 

The  best  correlations  of  experiment  and  prediction  are  for  relative 
thicker  laminates  impacted  by  7.5  mm  fragment  simulators  reported  by 
Simmons  et  al.  (10)  and  shown  as  relative  ballistic  limit  in  Fig.  7.  The  better 
correlation  may  be  related  to  the  higher  relative  target  thickness  to  projectile 
diameter  ratio  in  this  case,  which  ensures  that  the  first  stage  of  plug 
acceleration  dominates  the  problem  and  this  is  where  the  model  is  more 
accurate.  In  addition  with  the  smaller  diameter  projectile  plugging  is  less 
significant  and  this  is  the  aspect  where  the  model  poorly  represents  the 
simulation  of  behaviour.  This  type  of  problem,  thick  targets  and  small 
fragments,  is  the  more  usual  and  therefore  the  more  useful  application  in 
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practice.  In  addition  it  represents  the  class  of  problem  where  organic 
composites,  to  which  it  may  be  possible  to  extend  the  approach,  are  usefully 
applied. 


The  model  can  thus  be  used  with  metallic  laminates  to  estimate 
overall  ballistic  performance.  It  can  also  be  used  to  study  the  effects  of 
variations  in  laminate  structure  on  expected  performance,  parametric  studies, 
provided  one  is  aware  of  how  the  limiting  assumptions  may  influence  the 
output.  Thus  a  radical  change  in  material  type  from  a  ductile  to  very  brittle 
material  or  from  a  bonded  to  an  unbonded  laminate  may  result  in  changes  which 
the  model  does  not  include.  However,  number  of  layers,  bond  strength, 
laminate  thickness,  projectile  dimensions  and  some  variations  in  alloy  type  are 
all  parameters  amenable  to  study.  The  model  also  gives  a  picture  of  the 
physics  of  the  deformation  and  what  are  the  appropriate  properties  to  measure 
which  will  determine  ballistic  performance. 

The  computer  program,  LAMP,  is  listed  in  Appendix  1  with  a  typical 
input  and  output  as  well  as  a  discussion  of  the  output  parameters.  The 
program  was  written  to  estimate  ballistic  limits  by  incrementing  the  impact 
velocity  in  5  ms-1  steps.  The  understanding  of  behaviour  is  not  sufficient  to 
place  much  value  on  an  effort  to  predict  output  velocities  above  the  ballistic 
limit  using  such  a  model,  when  this  can  be  done  as  efficiently  and  effectively 
using  other  techniques  [91.  The  program  stops  and  prints  output  parameters  as 
soon  as  the  ballistic  limit  is  exceeded.  In  some  cases  the  residual  velocity 
predicted  can  be  substantial  as  the  small  incremental  velocity  steps  can  result 
in  a  sudden  change  in  one  layer  thickness  to  be  perforated  in  stage  n.  This  is 
as  would  be  expected  from  such  a  model  and  does  not  make  the  predicted 
ballistic  limit  of  any  less  value.  What  it  does  do  is  highlight  the  importance  of 
the  model  being  for  multi-layers  of  which  four  is  about  the  lower  limit.  A 
simple  two  layer  structure  can  be  handled  by  a  simple  procedure  developed 
earlier  [4],  By  combining  these  tools  one  has  a  useful  approach  to  studying  the 
ballistic  behaviour  of  laminates  which  should  reduce  considerably  the 
requirement  for  ballistic  testing.  In  addition  the  test  procedures  developed 
provide  simple  tools  for  studying  those  laminate  deformation  mechanisms  of 
importance  for  ballistic  performance. 


5.  CONCLUSION 

The  perforation  of  metallic  laminates  has  been  modelled  by  a  two  stage  process 
of  plug  acceleration  followed  by  dishing  behaviour.  Test  procedures  to 
determine  laminate  characteristics  were  outlined,  and  computations  compared 
with  empirical  data  for  ballistic  limits  for  a  range  of  laminates.  The  model 
makes  reasonable  predictions  of  ballistic  limit,  however,  its  main  use  is  shown 
to  be  as  a  tool  for  parametric  BtudieS  of  behaviour  in  relation  to  laminate 
design.  Limitations  in  the  application  of  the  model  relate  to  effects  of 
bonding  and  material  ductility  on  the  ductile/brittle  failure  response,  under 
impact  conditions. 
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Table  I  Metallic  Laminates  -  Physical  and  Mechanical  Properties 


PROPERTY 

L9 

LI  2 

L10 

MATERIAL 

L28 

CODE 

D1 

Ml 

M2 

M3 

Alloy 

7075 

7075 

7075 

7075 

2024 

2024 

7075 

2090 

No.  of  Layers 

6 

6 

6 

6 

4 

4 

6 

4 

Adhesive 

HYSOL** 

HYSOL 

HYSOL 

HYSOL 

- 

FM73 

FM73 

FM7  3 

Kevlar  Present? 

(Y/N) 

N 

Y 

N 

N 

N 

N 

N 

N 

Total  Thickness* (mm) 

6.38 

8.17 

7.98 

8.17 

6.36 

7.00 

7.40 

6.80 

Lamine  Thickness  (mm) 

1.02 

1.02 

1.02 

1.02 

1.59 

1.59 

1.02 

1.59 

Mean  Bond-line  Thickness* 
(lan) 

57 

414 

366 

414 

- 

213 

252 

147 

Vol  Fraction,  Alloy* 

0.96 

0.75 

0.77 

0.75 

1.0 

0.91 

0.83 

0.94 

Areal  Density^ 

(kg  m~2) 

17.14 

19.43 

18.90 

19.21 

17.81 

18.00 

18.58 

17.02 

Volume  Density^ 

(kg  m~3) 

2687 

2379 

2385 

2353 

2800 

2570 

2511 

2503 

0.2%  Proof  St^(HPa) 

483 

407 

387 

377 

355 

322 

418 

494 

uts'ihp") 

540 

511 

433 

422 

595 

473 

466 

517 

%  Elongation^ (% ) 

212 

-12 

212 

212 

-13 

213 

212 

29 

Note8i  *  Typical  values 

f  Calculated  values  based  upon  known  densities  of  constituents  and 
typical  values  (*) 

1  Predicted  values  based  on  Rules  of  Mixture  and  assuming  negligible 
tensile  strength  for  adhesive 
**  Hysol  EA9309.3  (NA) 


18 


Table  H 


Laminate  Input  Data  for  Use  in  Program  "LAMP' 


D1 

Ml 

MATERIAL  CODE 

L9  L12 

L28 

M2 

M3 

Thickness  (mm) 

6.36 

7.00 

6.42 

8.20 

8.50 

7.38 

6.80 

Density  (kg  m  ) 

2800 

2570 

2690 

2380 

2310 

2511 

2503 

(MPa) 

3020 

3715 

6920 

13335 

12302 

8913 

4753 

nl 

0.34 

0.43 

0.55 

0.82 

0.88 

0.66 

0.42 

eo 

.05 

.05 

.05 

.05 

.05 

.05 

.05 

o2  (MPa) 

3020 

2042 

2570 

2240 

2820 

2188 

1875 

n2 

0.34 

0.29 

0.21 

0.214 

0.38 

0.20 

0.11 

c  (MPa) 
z 

180.0 

188. 

178. 

166. 

134. 

164. 

195. 

op  (MPa) 

332. 

411. 

516. 

433. 

445. 

512. 

536. 

*c(*> 

0.1 

60 

2 

9 

50 

40 

32 

ou  (MPa) 

514 

412 

538 

503 

393 

487 

514 
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Table  in  Terminal  Ballistic  Data  for  a  Selection  of  Experimental 
Laminates 


Material 

Material  Type 

V-50  Values  (m/s) 

Code 

No  of 
plies 

Alloy 

Adhesive 

Experimental 

LAMP 

prediction 

D1 

4 

2024 

- 

-200 

235 

Ml 

4 

2024 

FM73 

-208 

285 

L9 

6 

7075 

Hysol 

210 

285 

LI  2 

6 

7075 

Hysol 
+  Kevlar 

236 

325 

L28 

6 

7075 

Hysol 

273 

305 

M2 

6 

7075 

FM73 

-220 

310 

M3 

4 

2090 

FM73 

189 

335 

Note:  Projectile  :  12.7  mm  diameter,  12.5  g  hard-steel  cylinder 
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APPENDIX 


The  program  for  computing  the  ballistic  limit  of  metallic  laminates 
is  called  LAMP.  It  takes  input  from  the  file  LIMP.DAT  and  output  is  provided 
in  a  file  LOMP.DAT.  The  program  increments  velocity  in  steps  of  5  ms-1  till  a 
velocity  is  reached  at  which  the  projectile  perforates  the  target,  defining  the 
ballistic  limit.  Table  A1  below  lists  a  typical  input  for  a  projectile  of  mass 
12.5  g  and  diameter  7.5  mm  and  using  material  L12.  The  output  is  given  at 
Table  A2  with  an  explanation  of  symbols  and  then  the  program  is  listed. 


Table  A1  T.mting  of  Input  Data  File  LIMP.DAT 


Parameter  Symbol  in  Text 

Symbol  in  Program 

Value 

Units 

Density 

P 

RO 

2380. 

kg  m~3 

Strength 

SOK1 

13335. 

MPa 

Exponent 

nl 

EX1 

0.82 

- 

Strain  Constant 

60 

EO 

0.05 

- 

Strength 

°2 

SOK2 

2240. 

- 

Exponent 

n2 

EX2 

0.214 

- 

Shear  Stress 

rs 

SS 

166. 

MPa 

Fibre  Stress 

°F 

BSG 

433. 

MPa 

Ultimate  Stress 

AB 

9. 

<°> 

Critical  Bend  Angle 

e 

c 

TS 

503. 

MPa 

Thickness 

ho 

HO 

.0082 

m 

No  of  Layers 

n 

NO 

6 

- 

Projectile  Mass 

m 

ASS 

.0127 

kg 

Projectile  Diameter 

D 

DIAM 

.0075 

m 

21 


Table  A2* 

Listing  of  Output 

File- 

LOMP.DAT 

BA  CL 

HT  STRN 

WK 

W 

VN 

V 

T 

.0013  .0069 

THROUGH 

.0075  0.089 

56.6 

87.3 

204.5 

203.7 

.00594 

W 

WP 

VR 

LNE 

LT 

LB 

357.7 

270.4 

41.6 

5 

2 

3 

EN 

VO=BALLISTIC  LIMIT 

350.7  235.0 

AREAL  DENSITY  KG/M2 
19.52 

TOTAL  THICKNESS  MM 

8.20 

NO.  LAYERS 
6 

LAYER  THICKNESS  MM 
1.367 


Symbols  :  (BA.CL.HT)  =  (Lg,  L^.h)  See  Figure  2. 

STRN,  Uniaxial  Strain  in  plug. 

WK,  Compressive  Work  Stage  I  (WN  in  text). 

W,  First  Line  -  Total  Work  Stage  I,  Second  Line  -  Total  Work  of 
Perforation. 

VN,  V  Plug  Velocity  at  end  of  Stage  I,  Projectile  Velocity  at  End  of 
Stage  I. 

T,  time  of  first  stage. 

VR,  residual  velocity  after  perforation. 

LNE,  Total  layers  to  be  perforated  in  Stage  n. 

(LT,LB)  =  (nj,  n2>  See  text. 

EN,  projectile  initial  kinetic  energy. 


LISTING  OF  PROGRAM  LAMP 


C  LAMP.FOR 

C  PROGRAM  TO  CALCULATE  PENETRATION  PARAMETERS  FOR 
MULTI-LAYER 

C  COMPOSITE  TARGETS  R.L. WOODWARD  AND  I.G.CROUCH.DSTO 
C  MATERIALS  RESEARCH  LABORATORY,  1988. 
OPEN(UNIT»l,FILE=’LIMP’)STATUS=’OLD’) 
OPEN(UNIT=2,FILE=’LOMP\STATUS=’NEW’) 
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21  FORMATYF5.0/F6.1/F5.3/F8.5/F6.1/F5.3/F6.1/F6.1/F5.2/F6.1 

1/F7.4/I4/F7.6/F7.6) 

READ(l,21)RO>SOKl,EXl,EO,SOK2,EX2,SS,BSG,AB,TS>HO,NO,ASS>DIAM 

75  FORM  AT(T4, ’BATll, ’CLT18,’HTT25,’STRNT33,’WKT44,’W’ 
2T55,’VNT64,’VT72,’TT80) 

WRITE(2,75) 

76  FORMAT(4X,F10.1,4XtF10.1,4X,F10.1,4X,I3)4X)I3>8X,I3) 

77  FORMAT(2X,F5.4,2X,F5.4,2X,F5.4,2X,F6.3,2X,F6.1,5X,F6.1, 
35X,F5.1,2X,F7.1,1X,F6.5) 

78  FORMAT(T3, THROUGH’) 

86  FORM  AT(T7 ,  ’  W  T22 ,’  W  GT38 ,  *VRT4  7 ,  ’LNE’) 

87  FORMAT<T7,’WT22,’WPT38,’VRT47,’LNET55,’LTT65,’LB’) 

88  FORMAT(T7,’EN’T20,'VO=BALLISTIC  LIMIT) 

79  FORMAT(2X,F10.1,10X,F6.1) 

51  FORMAT(T20,’AREAL  DENSITY  KG/ M2’) 

52  FORMAT(22X,F10.2) 

53  FORM AT(T20, TOT AL  THICHNESS  MM’) 

54  FORMAT(22X,F8.2) 

55  FORMAT(T20,’NO.  LAYERS’) 

56  FORMAT(22X,I3) 

57  FORM ATXT20, ’LAYER  THICKNESS  MM’) 

58  FORMAT(22X,F8.3) 

SOKl=SOKl  * 1000000 . 

SOK2=SOK2» 1000000 . 

SS=SS* 1000000. 

BSG=BSG»1000000. 

AB=.008727*AB 

TS=TS*1000000. 

TL=HO/NO 
ARED=RO  *  HO 
HOMM=HO*1000. 

TLMM=TL*1000. 

RAD=DIAM/2. 

Bl=((SOKl*EXl/RO)**.5)*2/(EXl+l) 

B2=((SOK2*EX2/RO)**.5)*2/(EX2+l) 

VOL=3.1416*RAD**2*HO 

VO=10. 

DVO=5. 

C  START  VEL.  INCREMENTING  IN  5  M/S  STEPS  TO  GET  BALL.  LIMIT 
DO  207  M=l,300 
EN=.5*ASS*(VO**2) 

ST=.001 

STRN=.001 

WK=0 

BA=0 

C  START  FIRST  STAGE  -  PLUG  ACCELERATION 
DO  102  L=l,5000 
IF(STRN.GE.EO)GO  TO  33 
VN=Bl*STRN**((EXl+l)/2) 
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F=VOL*SOKl  *(STRN  *  *EX1)  /HO 
GO  TO  34 

33  VN=Bl«EO**((EXl+l)/2)+B2*STRN**((EX2+l)/2)-B2*EO**((EX2+l>/2) 
F=VOL*SOK2*(STRN**EX2)/HO 

34  HT=HO*EXP(-STRN) 

T=<HO-HT)*2/VO 

FR=VN»T/2 

BAl=<VO+VN)*T/2 

DBA=BA1-BA 

BA=BA1 

CL=HO-BA 

WK=WK+F*DBA 

WS=3.14l59*RAD*SS*(HT»(HO-HT)+FR*(2*HT-FR» 

WF=(3.14159*SS*RAD*BA*»2) 

W=WS+WF+WK 
V2=VO**2-2*W/ASS 
IF(V2.LT.0.)GO  TO  206 
V=V2**.5 
T=T*1000 

IF(V.LE.VN)GO  TO  99 
IF(CL.LE.O)GO  TO  100 
STRN=STRN+ST 
102  CONTINUE 

C  BEGIN  SECOND  STAGE  -  DISHING 
99  LNI=1.+BA/TL 
TE=HO-<LNI*TL) 

LNE=NO-LNI 

SASS=ASS+3.1416*(RAD**2)*HT*RO 
RKE=0.5*SASS»(VN*  «2) 

WT=.3927*(DIAM**2)*HT*TS 
IF(AB.LT..004)GO  TO  308 
AB=AB*((HO/TE)**2) 

IF( AB.LT. 1 .5708)GO  TO  306 
AB=1.5708 

C  WORK-GOOD  BONDING 
306  WBG=.7863*DIAM*(TE**2)*BSG*AB 

LT=LNE/2 
LTD=2*LT 

IF(LTD.EQ.LNE)GO  TO  204 

LB=LT+1 

TT=TL*LT 

TB=TL*LB 

GO  TO  205 

204  LB=LT 
TT=LT*TL 
TB=TT 

C  WORK-LAMINATE  BOND  BROKEN 

205  WBP=.7853*DIAM*BSG*«TT**2)+(TB**2»»(1.5708-AB) 

GO  TO  320 
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C  WORK-NO  INTERLAMINAR  BONDING 
308  WBN=1.2335*DIAM*BSG*(TE**2)/LNE 
320  W2=W+WT+WBG+WBP+WBN 
WP=WT+WBG+WBP+WBN 
V2=VN  *  *2-2*  WP/SASS 
IF(V2.LT.0.)GO  TO  206 
VR=SQRT(V2) 

WRITE(2,77)BA,CL,HT,STRN,WK,W>VN,V',T 

WRITE<2,78) 

WRITE(2,87) 

WRITE(2,76)W2,WP,VR,LNE,LT,LB 
GO  TO  350 

206  VO=VO+DVO 

207  CONTINUE 

100  WRITE(2,101) 

101  FORMAT(T3,’CL=0  OR  OUT  OF  TIME’) 

350  WRITE(2,88) 

WRITE(2,79)EN,VO 
WRITE*  2, 51) 

WRITE(2,52)ARED 

WRITE<2,53) 

WRITE(2,54)HOMM 

WRTTE(2,55) 

WRITE(2,56)NO 

WRITE(2,57) 

WRITE(2,58)TLMM 

STOP 

END 
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(a) 


(b) 


(c) 


Figure  1  Schematic  of  the  perforation  process  in  a  laminated  target. 

(a)  The  initial  impact  conditions, 

(b)  Stage  I  where  the  projectile  moves  faster  than  the  target  plate 
material  in  front,  so  that  compressive  deformation  of  the  target 
material  occurs  as  well  as  shearing  between  the  material  ahead  of 
the  projectile  and  the  rest  of  the  target. 

(c)  Stage  U  when  the  projectile  and  * 'plug "  are  moving  at  the  same 
velocity,  delamination  of  unsheared  laminae  at  the  back  of  the 
target  occurs  and  deceleration  is  affected  by  resistance  of  the 
dishing  target  material. 
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fa) 


(b) 


Figure  2  (a)  The  deformation  of  a  laminated  target  at  the  end  of  Stage  I,  where 

LB  defines  the  position  of  the  back  of  the  plug  or  front  of  the 
projectile.  The  material  to  be  dished  in  Stage  E  is  determined  by 
the  next  interface  below  the  projectile  position,  thickness  TB> 
Other  symbols  relate  the  deformation  to  that  of  a  homogeneous 
target  (11. 

(b)  Deformation  of  the  laminate  in  Stage  E showing  the  debonding  and 
those  laminae  which  have  failed  by  shearing. 
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Figure  3  (a)  Typical  uniaxial  stress/uniaxial  strain  data  for  6-ply  7075  aluminium 

alloy  laminates  of  varying  adhesive  thickness  from  constrained 
compression  tests. 

(b)  Sample  of  a  6-ply  7075  aluminium  alloy  laminate  plate  which  has 

been  compressed  between  flat  dies  in  a  constrained  compression  test 
showing  the  typical  strain  distribution  and  the  radial  extrusion  of 
material  against  the  constraint  provided  by  the  surrounding  plate. 
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Figure  4  Load/displacement  test  data  for  monolithic  and  laminated  samples 

showing  the  construction  of  a  line  to  define  the  shear  strength  of  the 
laminate  using  the  area  under  the  load  displacement  curve  for  the 
laminate  and  its  through  thickness  as  fixed  parameters. 
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TOTAL  BEND  ANGLE.  Bldegrees) 


Figure  5  Typical  moment/angle  of  bend  data  for  laminates: 

A  -  unbonded,  B  &  C  -  poorly  bonded  and  D  -  well  bonded  -  the 
interlamellar  shear  strengths  measured  are  indicated  and  failure  in  shear 
by  a  sudden  drop  in  bending  moment  at  the  indicated  angles. 
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RELATIVE  V50  BALLISTIC  l 


Figure  7 


NUMBER  OF  LAYERS 


Comparison  of  model  calculations  and  experiment  (10)  for  7.5  mm 
fragment  simulating  projectiles  impacting  laminates  of  varying  thickness. 


SYMBOLS 


Longitudinal  wave  speed 
Projectile  diameter 
Force 

Laminate  initial  thickness,  plug  thickness  after  acceleration 

Displacement  of  back  (impact  side),  and  front  of  plug 
respectively 

Bending  moment  measured  in  four  point  bend  test 
Projectile  mass 

Number  of  layers  in  laminate,  number  of  layers  to  be 
penetrated  in  Stage  n 

Numbers  of  layers  in  two  sections  following  bond  failure  in 
Stage  II 

Thickness  of  laminate  to  be  perforated  in  Stage  II 

Thickness  of  one  layer  of  laminate 

Time 

Velocity  of  plug  after  acceleration,  of  projectile  after 
deceleration,  and  impact  velocity  respectively 

Velocity 

Work,  Work  done  in  plug  compression,  Work  of  shear,  Work  of 
WB.  W rji  friction,  Work  of  bending  and  Work  of  tensile 
deformations  respectively 

Width  of  sample  in  four  point  bend  test 

Yield  stress 

Stress,  ultimate  stress,  fibre  stress  respectively 

Curve  fitting  parameters,  including  combined  effect  of  flow 
stress  and  constraint,  equations  (1) 


Strain  and  strain  constant  curve  fitting  parameter  for 
equations  (1) 

Density 

Shear  flow  stress 

Angle  of  bend,  critical  angle  of  bend  for  bend  failure  in  four 
point  bend  test 

Curve  fitting  parameters  for  work  hardening  exponent, 
equation  (1) 
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